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Abstract 

We consider a gradient estimate for a conductivity problem whose inclusions are 
two neighboring insulators in three dimensions. When inclusions with an extreme 
conductivity (insulators or perfect conductors) are closely located, the gradient can be 
concentrated in between inclusions and then becomes arbitrarily large as the distance 
between inclusions approaches zero. The gradient estimate in between insulators in 
three dimensions has been regarded as a challenging problem, while the optimal blow¬ 
up rates in terms of the distance were successfully obtained for the other extreme 
conductivity problems in two and three dimensions, and are attained on the shortest 
line segment between inclusions. In this paper, we establish upper and lower bounds 
of gradients on the shortest line segment between two insulating unit spheres in three 
dimensions. These bounds present the optimal blow-up rate of gradient on the line 
segment which is substantially different from the rates in the other problems. 

1 Introduction 

Let Bi and B2 be bounded and simply connected domains in d = 2,3. We consider 
the following conductivity problem: for a given harmonic function H defined in 

fv • ((fc-1 )x(b^uB2) + 1 ) = 0 

'|u(x) — H{x) = as |x| —>• 00, 

where y is the characteristic function. Two inclusions Bi and B2 are conductors with 
conductivity k ^ 1, embedded in the background with conductivity 1. For an extreme 
conductivity A: = 0 or 00, a modified model has been used, see GH) and m- Let e 
denote the distance between Bi and B2, i.e., 

e := dist(Bi, i?2), 

and we assume that the distance e is small. 

The problem is to estimate |Vm| in the narrow region in between inclusions. This 
was raised by Babuska in relation to the study of material failure of composites [4] . In 
fiber-reinforced composites which consist of stiff fibrous inclusions and the matrix, high 
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shear stress concentrations can occur in between closely spaced neighboring inclusions. 
It is important to estimate the shear stress tensor Vu, while u means the out-of-plane 
displacement, and the inclusions Bi and B2 are the cross-sections of fibers. Many 
studies have been extensively conducted on the gradient estimate due to such practical 
importance. 

Successful results have been achieved in all cases except three dimensional insula¬ 
tors which we consider in this paper. Such successful results can be divided into three 
cases when k stays away from 0 and 00, when k degenerates to either 0 (insulating) 
or 00 (perfectly conducting) in two dimensions, and when fc = 00 in three and higher 
dimensions. On the other hand, this paper deals with the exceptional case when k = 0 
in three dimensions. We prove the occurrence of concentration in the narrow region, 
and also established the optimal blow-up rate for |Vu| on the shortest line segment 
between two insulating unit spheres in terms of e. 

We give a brief description of three successful cases mentioned above. In the first 
case when k stays away from 0 and 00, i.e., ci < k < C2 for some positive constants 
Cl and C2, it was proved by Li-Vogelius m that |Vu| remains bounded regardless of 
the distance e between inclusions. The boundedness result was extended to elliptic 
systems by Li-Nirenberg [ 15 ]. 

In the second case when k is either 0 or 00 in the two dimensional problem, the 
gradient Vu can become unbounded as the distance e tends to 0 , and the generic 
blow-up rate of |Vm| is l/Ve For two circular inclusions, the blow-up rate l/\/e was 
derived by Budiansky-Carrier [ 7 ], and Kang-Lim et.al. Elia established the precise 
dependence of |Vu| on e, radii of disks and k € [ 0 ,oo]. For inclusions in a sufficiently 
general class of shapes in it was shown by Yun [ 20 l[ 21 ] that the blow-up rate l/\/e 
is valid at fc = 0 or 00, see also [ 18 ] for an enhancement of concentration. Taking it a 
step further, an asymptotic for the distribution Vrt was established by Kang-Lim-Yun 
E), when Bi and B2 are disks. For sufficiently general shapes of inclusions, recent 
results by Ammari et.al. [I] and Kang-Lee-Yun [ 5 ] yield a numerically stable method 
to well describe the asymptotic behavior of Vu in R^. 

In the third case when /c = 00 in three and higher dimensions, Bao-Li-Yin [^ proved 
that the generic blow-up rate for the perfectly conducting inclusions is |eloge|“^ in 
three dimensions and |e|~^ in higher ones, see also | 5 | for multiple inclusions. Lim-Yun 
m also found the explicit dependency of |Vu| on the radii as well as the distance e, 
when two inclusions are spheres in three and higher dimensions, see also [nilllKla]. 
Morever, an asymptotic for the distribution Vu was established by Kang-Lim-Yun 
nni and Lim-Yu m when Bi and B2 are two perfectly conducting spheres in R^. 

It is important to estimate |Vu| on the shortest line segment between inclusions. 
In the second and third cases mentioned above, the asymptotic behaviors of Vm show 
in 0 [iiiHiiiniiis] that the generic blow-up rates can be attained on the shortest 
line segment. Since high concentration of Vw results from the reflections repeated 
infinitely between inclusions, refer to nnnKni, it is reasonable to estimate |Vu| on 
the shortest line segment to obtain the blow-up rate. 

This paper is mainly concerned with the exceptional case when the inclusions is 
insulators (k = 0 ) in three dimensions. The gradient estimate in this case has been 
regarded as a challenging problem. An upper bound of |Vu| with order l/i/e was 
derived by Bao-Li-Yin [B] in this problem. To our best knowledge, there has not been 
any updated or improved result yet. 

In this paper, we prove that the optimal blow-up rate of |Vrt| is I/e 2 on the 
shortest line segment between two insulating unit spheres in three dimensions. To do 
so, we establish the upper and lower bounds of |Vu| on the line segment. It is worthy 
mentioning that the blow-up rate is substantially different from the rates known in 
the other extreme conductivity problems. In terms of methodology, a new method is 
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employed to derive the optimal blow-up rate due to the different nature of the problem. 


2 Main Results 

In this paper, we assume that Bi and B 2 are a pair of unit spheres e apart as follows: 

Si =Bi (- 1-1 0 , 0 ) andS 2 =Si (1 + 1 , 0 , 0 ) . 

The quantity e means the distance between Bi and S 2 , and is supposed to be small 
as mentioned in the introduction. The centers of Bi and B 2 lie on the x axis. For any 
harmonic function H defined on let u be the solution to the following conductivity 
problem whose conductivity k degenerates to 0: 

{ Au = 0 in \ (Si U S 2 ), 

d^u = 0 on d{Bi U S 2 ), 

m(x)-S( x) = O (^ 1 ^^ as|x|-^oo. 

Here, x = (x, y, z) in R^. 

This paper has two main results that provide an optimal blow-up rate 1/e 2 for 
|Vu| on the shortest line segment between two insulating Si and S 2 . 

Theorem 2.1 (Upper Bound) For any harmonic function FI defined in R^, let u 
he the solution to (EH). There is a constant C independent of the distance e > 0 such 
that 

|Vu(x,0,0)|<U^ 
e 2 

for |x| < I and any small e > 0. 

The proof of this theorem is presented immediately after Proposition 12.51 since it is 
based on Propositions 12.3112.41 and 12.51 

Theorem 2.2 (Lower Bound) Assume that 

H{x,y,z) = y 

in R^. Let u be the solution to EH for S. There is a constant C independent of the 
distance e > 0 such that 

|Vn(x,0,0)|>C^ 
e 2 

for |x| < I and any small e > 0. 

The proof of Theorem 12.21 shall be provided in Section [S] by modifying the proof of 
Proposition 12.51 

To show Theorem EH we use Propositions 12.3112.41 and 12.51 The first proposition 
means that the directional derivative dxU is bounded regardless of e, where the x axis 
passes through two centers of Hi and H 2 . 

Proposition 2.3 Let u be the solution to EH for a harmonic function F[ in R^ as 
given in Theorem \2.1[ There is a constant C independent of e > 0 such that 

|dxM(x, 0,0)1 < C 

for \x\ < I and any small e > 0. 
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This proof of Proposition 12.31 is provided in Subsection EH 

To estimate the other directional derivatives, we use an upper bound in the second 
proposition which simplifies this problem. 

Proposition 2.4 For a harmonic function H defined in let u he the solutions to 
for H. Then, there exists a large constant M independent of e > 0 such that for 
|x| < I and e > 0, 

0,0)1 < |9y'Ui(a::, 0,0)| , 

where ui is the solutions to (EB for Hi{x,y,z) = My. 

In Subsection 13.21 we prove Proposition 12.41 

The following proposition is an essential part of this paper that actually yields the 
first main result. 

Proposition 2.5 Let u be the solution to (EH) for H[x,y,z) = y. Then, there is a 
constant C independent of e > 0 such that 

\Vu{x,0,0)\<C^ 
e 2 

for |a;| < | and any small e > 0. 

In Section H we prove Proposition 12.51 To do so, we present Proposition 14.11 that 
obviously implies Proposition 12.51 Thus, Section |4] is mainly devoted to showing 
Proposition 14.11 

Now, we are ready to prove the first main result by virtue of three propositions 
above. 

Proof of Theorem 12.II . Propositions 12.31 implies that the directional derivative dxU 
is bounded independently of e > 0 on the line segment between Bi and i? 2 , supposed 
that the distance e is small enough. Next, we consider the other directional derivatives. 
By Propositions 12.41 and 12.51 

|ayu(a:, 0,0)1 < |VMi(a:, 0,0)1 < Ci 

e 2 

for |a;| < | and any small e > 0, where ui is the solution to (12.ip for H = Miy for a 
large Mi > 0. Similarly, we can choose a positive constant M 2 independent of e > 0 
so that 

0,0)1 < |Vu2(x,0,0)1 < C2^^ 

e 2 

for \x\ < I and any small e > 0 where U 2 is the solution to (12.11) for H = M 2 Z for a 
large M 2 > 0. Therefore, we complete the proof. □ 

3 Representations of the solution u 

Two representations m and (13.5p of the solution u are introduced in this section. 
The first representation is used to prove Proposition 12.51 and the second one is for 
Propositions 12.31 and 12.41 whose proofs are presented in Subsections 13.11 and 13.21 re¬ 
spectively. 

First, the solution can be decomposed into three harmonic functions as 

u{x,y,z) = H{x,y,z) + RBi[x,y,z) F RbAx.Vm), (3.1) 
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where the harmonic function Rb^ is defined uniquely in \ Bi and satisfies the decay 
condition = O as |x| —)> oo for i = 1 , 2. The decomposition can be derived 

from the representation of u as a sum of H and two single layer potentials on dBi 
and on 552, respectively. For details, refer to the invertibility of —\l — K* presented 
in Section 2 in [1], where K* is the the Neumann-Ponicare operator. It is worthy 
mentioning that Rbi is the reflection oi H + Rb^ occurring on the insulating inclusion 
Bi only, and similarly Rb^ is the reflection oi H -\- Rb^ on B 2 . 

In this paper two harmonic functions Rb^ and Rb 2 play a important rule. They 
are used for proving ProDOsition l2.5l that is actually the first main result in this paper. 
We study the properties of Rb^ and Rb^ in Section S] where the proof of Proposition 
12.51 is presented. 

Second, another representation of u is also introduced in Lemma l3.II This involves 
the derivations of Propositions 12.31 and 12.41 To illustrate the second representation, 
we consider the reflection only for a single inclusion Bq that denotes the unit sphere 
with center (0,0,0), i.e., 

Bo = Si (0,0,0). 

For any harmonic function h defined in a neighborhood containing Bq, let Ro{h) be 
the reflection of h with respect to Bq, i.e., 

{ ARo{h)=0 

di,{h + Ro{h)) = 0 ondBo, 

i?o(h)(x) = O as |x| -)> 00. 

In the spherical coordinate system, 

Ro{h){p,0,4>) = - J h{s,0,4>)ds (3.3) 

for p > 1, where {x,y,z) = {p cos 0 sin cj), p sin 0 sin cj), p cos cj)) in In the Cartesian 
coordinate system, 

— ,0, Oj—— j sdyh{s,0,0)ds (3.4) 

X J X Jq 

for X >1. 

Similarly, for any harmonic function h defined in a neighborhood containing Bi, 
we define Ri{h) as the reflection of h with respect to for i = 1, 2 as follows: 


in R.^ \ Bi, 
on dBi, 

as |x| —>• oo. 

Lemma 3.1 The following sum converges to the solution u in the sense of the Sovolev 
space IF^’°°(R^ \ (Bi U B 2 )), i.e., 

u(x) =B(x) + Bi(B)(x) + B 2 (B)(x) 

00 

+ ^(BiB2)”(B)(x) + (BiB2)"Bi(B)(x) 

n—1 
00 

+ ^(B2Bi)"(B)(x) + (B2Bi)"B2(B)(x) (3.5) 

n—1 


AR,(h) = 0 
di,{h + Riih)) = 0 

Ri{h){x) = O 


for any x G R^ \ (Bi U B 2 ). 


5 




Proof. We begin by showing two properties (EH) and (HTfl) that are essential to 
prove the convergence of the series (1^ . Let h he a harmonic function defined in 
(K^ \ Bi) with the decay condition /i(x) = 0( j^) as |x| —>■ oo. 

We first show that 



R2{h)d, {R2{h)) dS 


1 

- (TT^ 


hd^h dS 


IdBi 


(3.6) 


Since Bi is the unit sphere with the center (—1 — |,0, O), the function h can be 
expressed in terms of spherical harmonic functions whose center is (—1 — §, 0, O). By 
the decay condition of h, we have 


[ hd^h dS 

> 

f 

hd^h dS 

JdBi 


d|x-(-l-ie,0,0)| = l+£ 



= 

/ 

iV/rp 

dV 



7|x-(-l-ie,0,0)|>l+e 



> 

/ hd^h dS 





JdBi 




> 

/ R2{h){d,R2{h)) dS 




JdBi 




Thus, we have (EH). 

Second, we prove the existence of a constant Ci such that 


Cl 


hd^h dS 


IdBi 


> \\R 2 m W^'“(R\B2) 


(3.7) 


The mean value property for harmonic functions yields the inequality 


f hdjih dS 

1 

2 

f |V/i|2 dV 

JdBi 


dR3\Bi 


— ^2 II V/l||^oo(|x_(i+ie,0,0)|<l+ie) 


>Cz{ max (^(x)) 

y|x-(l+ie,0,0)|<l + ie 


min 

|x-(l + ie.0,0)|<l+i 



The positive constants C 2 and are used above regardless of choosing a harmonic 
function h in \ Bi)- Note that R 2 {h) can be extended as a harmonic function 

into the area |x — (1 + ie,0, 0)| > ^1 , and that the analogous formula for i? 2 (/i) 

with (13.311 is valid in the extended domain. Thus, the formula for i? 2 (^) implies 


C 4 


hd^h dS 


laBi 


> max |i? 2 (/i)(x)|. 

|x-(l + ie.0,0)|>^^ 


The distance between the boundaries of B 2 and the extended domain of i? 2 (^) is 
greater than je. For any point in \ 7^2, the harmonic function R 2 {h) is defined in 
the open sphere whose center is the given point and the radius is je. Thus, a gradient 
estimate for harmonic functions implies 


C 5 



hd^h dS 


> \\^R2{h)h °°(R3\B2) ’ 
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and moreover, the bound (EH) for the higher order derivatives can also be derived in 
the same way. 

Now, we are ready to prove this lemma. Applying (lABl) repeatedly, we have 


'dBi 


{RiR2T{H)dA{RiR2T{H))dS 


IdBi 


{RiR2T{RiH)d, {{RiR2T{RiH)) dS 


idBz 


{R2Rir{H)d,{{R2RiT{H))dS 


'dB2 


{R2RiT{R2H)d, {{R2 RiT{R 2H)) dS 


< 


< 


_f_y 

(1 + £)3(2n-l) 

'' ' i—\ 


IdBi 


Ri{H)d,Ri{H)dS 


Hd^HdS 


(l + ,)3(2n-l)^ JOB. 

for any n = 1, 2, 3, • • •. By (1X71) . we have 

\\{R,R2nH)\\ 

+ \\{R2RinH)\\ + \\{R2Rir{R2H)\\ 

M/4.“(R3\B2) 

2 


- (1 + e)3(2n-l) X! 

^ ' i—1 


Hd^HdS 


IdBi 


(3.8) 


for any n = 1,2,3,---, where Ci is the constant in (EZl). This implies that the series 
in the right hand side of (1X51) are convergent in the sense of \ {Bi U B 2 )) 

and thus satisfies m- Hence, the series converges to the solution u. □ 


ProDOsitions l2.3l and l2.4l can be derived from basic properties of the representations 
introduced in this section. 


3.1 Proof of Proposition 12.31 


We begin in considering the case of a single inclusion Bq. As defined before, Ro{h) 
is the reflection of a given harmonic function h with respect to Bq. By the equation 

(IXX) . 

0,0) = —f i, 0,0 I for x > 1. (3.9) 

x'^ \x ) 

Here, ^ in the equation means the first coordinate of (y,0, O) at which the image 
charge of (x, 0, 0) with respect to Bq is located. 

Lemma IXTI implies that the solution u results from the recursive reflections on Bi 
and H 2 . Dealing with the recursive reflections, we define ri(x) and r 2 (x) as the first 
coordinates of the image charges of (x, 0,0) with respect to Bi and i? 2 , respectively. 
Thus, 


ri(x) = 


1 + i 


(- 1 ) 




for X > 
for X < -. 
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For |x| < |, we define two sequaneces (r^„) and (rBn) as 


rA2n-i = (?’ 2 ?'i)” r 2 ix) and rA2n = -(rir 2 )"(a:), 


and 


rB 2 n-i = -irir 2 )^~^ri{x) and rB 2 n = (?' 2 ?'i)"(a;) 
for n = 1,2,3,-where {rirj){x) = ri{rj{x)) for {i,j} = {1,2}. Applying (13.91) to 
(13.5p . 

dxu{x, 0,0) =dc,H{x, 0,0) 

oo / n \ ^ 

+ n (l + I - {{-ir+^TAn, 0, 0) 

n—1 \fc—1 / 

oo / n \ ^ 

+ E(-l)” n (l + I - d^H {{-irVBn, 0, 0) 

n—1 Vfc—1 / 

for any |x| < |. 

Indeed, two positive sequences (r^„) and {tbu) are increasing and converge to a 
number that is i/e + 0(e). There are some properties that can be shown easily: 

I < rAn{x) < 2-/e, I < rBn{x) < 2Ve 

for \x\ < I and n = 0,1, 2, 3, • • •, and 

^\/e < (r 2 ri)”(a:) and <-(rir 2 )"(a:) 

for any n > — §■ These properties follow immediately from Lemma |4.8I 

in this paper. 

Now, we prove the boundedness in Proposition [231 Let a a = — 1 and as = 1- For 
j = A or B, 

oo / n \ ^ 

E(-l)” n + 5,if(a,(-irr,„,0,0) 


n—1 \fc=l 


oo /2n—1 


< 


, 1 + 2 “ ’’Jfe 




< 


E n( 

n—1 \ k—1 

X \da:H{-ajrj(^ 2 h-i}, 0 , 0 ) - d^H{ajrj(^ 2 h), 0 , 0 )\ + { 6 y/e + O{e))\d^H{rj(^ 2 h), 0 , 0 ) 

E ('^V^II^^-^(^’®’*^)llL“([-2v^,2Vd) + (6\/e + 0(e)) ||9a:7J(a;, 0,0)||j;^oo([_2^,2Vd) 


[_ 20-\/e 


h—1 


X (4Ve ||cl^iL(a:, 0,0)||^„„j.j_2^ 2Vi]) (6v^ + ^(^)) (a::, 0,0)||^oo([_2^ 2^/^])) 

C (||^a;-^(2^^*^>*^)||L“([-2Ve,2^/£]) III1 l“([- 2 Ve,2Ve])) ’ 

Hence, we have the boundedness of |cla;M(a:, 0, 0)| as 

\\dxu{x, 0,0)||^oo([_ e 6]) < C (^\\d^H{x, 0,0) ||sooq_2^_2Vd) ^')\\l°°{[-2^,2^)^ 


< 


for small e > 0. 


□ 





3.2 Proof of Proposition 12.41 

In the following lemma, we first consider the model of a single inclusion Bq that is 
much simpler than our model of two inclusions. Second, applying the lemma to Lemma 
o we prove Proposition [mi 

Lemma 3.2 Let h be a harmonic function defined in a neighborhood containing the 
closure of the unit sphere Bq = i3i(0, 0,0), and let Ro{h) be defined as (13.21) . Suppose 
that 

dyh(x,0,0)>0 and 9a;9y/i(x, 0,0) > 0 

for 0 < a; < 1. Then, 

dyRo{h){x,0,0) > 0 and dxdyRo{h){x,0,0) < 0 

for a: > 1. 

Proof. By (13.41) and the assumption, 

0, 0) = , 0, 0 I — — f sdyh (s,0,0) ds 

x-^ \x J X Jo 

= — y s (^dyh , 0,0^ — dyh (s, 0, 0)^ ds > 0 


for a; > 1. Thus, we have the first bound. 

We can get the second bound from dsn as follows: 


x"^ ^ \x 


1 


1 


1 


dxdyRo{h){x,0,0) =-—dyh[ -,0,0] +-^ J sdyh{s,0,0) ds - —d^dyh [ -,0,0 


,.5 


1 


< — 


1 


s I Adyh I —, 0, 0 ) — dyh (s, 0, 0) ) ds < 0. 


□ 


Now, we are ready to prove Proposition [531 
Proof of Proposition 12.41 From the definitions of Ri and R 2 , 

dyRi{y){x, 0,0) = —-^ > 0, (3.10) 

{-iydxdyRfiy){x,Q,f)) = — - \ . > 0 (3.11) 

2(l + f-(-l)*a:)" 

for — 1 — I < (—l)*a: < | and * = 1, 2. There is a large M > 0 such that 
dyRfiH){x,0,0) < dyRfiMy){x,0,0) 

and 

{-iydxdyRfiH){x,0,0) < {-iydxdyRfiMy){x,0,0) 

for —1 — I < {—ifix < I and i = 1,2, since RfiMy) = MRfiy). By mathematical 
induction. Lemma 151^ allows the upper and lower bounds of dy{RiR2)^{H){x,0,0) 
and dxdy{RiR 2 )^{H){x, 0,0) and so on such that for any n = 1, 2, 3, 4, • • • , 

dy{RiR2T{H){x,Q,t)) < dy{RiR2){My){x,t),Q), 
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d,dy{RiR2nH){x,0,0) > d,dy{RiR2nMy){x,0,0) 

and 

dyiRiR2r-^RiiH)ix,0,0) < dy{RiR2r-^RiiMy){x,0,0), 
d,dy{RiR2r-^RiiH){x,0,0) > d,dyiRiR2r-^Ri{My){x,0,0) 
for—|<a;<l + |, and 

9y(i?2i?i)"(i?)(x,0,0) < dyiR2Rir{My){x,0,0), 
d,dyiR2Rir{H)ix,0,0) < d,dy{R2RinMy)ix,0,0) 

and 

dy{R2Rir-^RiiH){x,0,0) < dy{R2Rir-^R2iMy){x,0,0), 
d,dy{R2Rir-^RiiH){x,0,0) < d,dy{R2Rir-^R2{MyKx,0,0) 
for — 1 — I < a: < |. By Lemma [3?T1 we have the upper bound of dyu{x, 0,0) as 
dyu(x, 0,0) <M + dyRi{My){x, 0,0) + dyR 2 {My){x, 0,0) 

OO 

+ ^ dyiRiR2riMy)ix, 0, 0) + 0,0) 

n—1 

OO 

+ ^ dyiR2Rir{My){x, 0, 0) + dy{R2RiTR2{My){x, 0,0) 

n—1 

=dyUi{x, 0,0) 

for |x| < when m is the solutions to (EH) for Hi{x,y,z) = My. The lower bound 
is also obtained in the same way. 

□ 


4 Proof of Proposition 12.5 

In this section, we assume that 

Hix,y,z) = y inM.^ 

and u is the solution to (EH) for H = y. As defined in the decomposition (EH , two 
harmonic functions Rbi and Rb 2 satisfy 

u{x, y,z) = y + Rbi {x, y, z) + Rb^ {x, y, z), 

where i?_Bi is defined in \ Bi and satisfies the decay condition Rb^ = O as 

|x| —>• OO for f = 1, 2. 

This section is mainly devoted to proving Proposition 14.11 that obviously implies 
Proposition 12.51 since dzu{x,0,0) = 0 and dxu{x,0,0) is bounded for \x\ < | by 
Proposition 12.31 

Proposition 4.1 dyRBi (x, 0,0) is decreasing in [—|,oo) and 

o<a,i?B, (-|,o,o) <^- ( 4 . 1 ) 

6 2 

R follows from the property dyRBi (x, 0,0) = dyRB 2 (—x, 0,0) that 

dyu {x, 0,0) = 1 + dyRBi 0,0) + dyRB 2 {x, 0,0) < (4.2) 

e~ 

for |a;| < f. 
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Here and throughout this paper, oi < bi means oi < Cibi and 02 — 62 stands for 
^02 < ^2 < C 2 a 2 for some constants Ci and C 2 independent of e. 

The proof of Proposition 14.11 is presented in Subsection 14.21 based on the lemmas 
in Subsection KT[ 

4.1 Basic Properties of dyRBi{x,0,0) 

We consider the behavior of (x, 0,0) to derive Proposition 14.11 In this subsec¬ 

tion, H{x,y,z) = 2 / in as assumed early in Sectional For convenience, we define 
the function P : [ 1 , 00 ) —)► R as 

P{x) = dyRsi - 1 - 0, 0 ) 

that is a horizontal shift of dyRBi{x,0,0). The translation moves the left inclusion 
Bi to Pi(0,0,0) so that the domain of P is the interval [l,oo) with the initial point 
1. The symmetry between (x, 0, 0) and SyPsj (x, 0,0) yields 

dyu{x, 0,0) = dyH{x, 0,0) -I- dyRsi (x, 0,0) -I- dyRB 2 {x, 0,0) 

= 1-I-P^x-1-1-I-|^ -l-P ^-x + 1 + 1) 

for |x| < since H{x,y,z) = H{—x,y,z). We study the behavior of P(x) especially 
for small x — 1 > 0 to prove the bound 

Pix) < (4.3) 

e 2 

that means Propositions 14.II and 12.51 

The basic properties of P are introduced in this subsection. First, Lemma 14.21 
provides a fundamental equation (14.6p of P that yields almost all properties of P 
including the main result in this paper. Second, Lemmas 14.31 and 14.41 describe the 
geometric behavior of P. Third, Lemma 14.51 presents an estimate for the integral 
value of P that determines the blow-up rate of P as e approachs 0. Finally, based 
on these properties, our main result (14.31) can be obtained in Subsection 14.21 to show 
Proposition 14.11 

Dealing with the lemmas, it is necessary to define two special points pi and p 2 as 
a pair of solutions to 

2 -I- e - - = X, (4.4) 

X 

where pi < P 2 - Indeed, they are the fixed points of the composition of two Kelvin 
transforms with respect to Pi(0,0, 0) and Pi(2-|-e, 0, 0), and can be calculated directly 
as 

Pi = 1 + (—l)*\/e -|- Oi{e) (4-5) 

for f = 1 , 2 . 

In the following lemma, we establish a fundamental equation of P which is essential 
in deriving many properties of P in this paper. 

Lemma 4.2 The function P : [1, 00 ) —>■ R. satisfies 

( _! 

2 +sP{2 + e- s)ds + ^--^ = P{x) (4.6) 

X J X Jq 2 x'^ 

for any x > 1 . 
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Proof. We first recall the property (El) of the reflection with respect to a single 
inclusion. Second, we apply dH to the case of two neighboring inclusions thus to 
derive dUl). 

We first consider the reflection only for a single inclusion Bq that is the unit sphere 
i?i(0,0, 0). For any harmonic function h defined in a neighborhood containing i?o, the 
reflection Ro{h) with respect to Bq satisfies 

— ,0, Oj—— / sdyh{s,0,0)ds (4.7) 

X ) X Jq 

for a; > 1, that is actually (El- 

Second, we consider the solution u to (12.11) . It can be decomposed into three 
harmonic functions as u{x, y, z) = H{x, y, z) + Rb^ (x, y, z) + Rb^ (x, y, z). From def¬ 
inition, Rbi {x, y, z) can be regarded as the reflection oi h = H (x, y, z) + Rb 2 (x, y, z) 
with respect to i?i = So + (—1 — f, 0, O), and RB 2 {x,y, z) = RBi{—x,y,z) due to 
the symmetric property of H{x,y, z) = y. Since P{x) = dyRBi {x — 1 — |, 0, O), the 
equality (El thus yields 


P{x)=dyRB, (a;-l-|, 0 , 0 ) 


= dyRo 


(-( 


x-1- -,y,z 


^ +Rb2{x-1- (a::, 0,0) 


= dyRt 


h (y + Rbi (^-x + 1 + ^,y,z'j'j (x,0,0) 


1 1 


- ^ + ^9yRB, Tx + ^ + 2’j ■ ^ ■ X 


1 1 


sd, 


iRbi ( 


— 5 - 1-1 + —, 0 , 0 ) (is 


1 1 1 

2 x^ x^ 


P 2 + e- 


sP{2 + e — s)ds. 


Thus, we have this lemma. 


□ 


Lemma describes the graph of P{x) as an application of Lemma [4.21 Thus, P{x) 
and dyRBi (x — 1 — 0, O) is positive and decreasing for x > 1. 

Lemma 4.3 The function P{x) satisfies 

P{x) > 0, P'(x) < 0, P"(x) > 0, P"'(x) < 0, P""(x) > 0 


for X > 1 and 

lim P(x) = 0. 

X—^OO 

Proof To prove the decay of P at infinity, we use the limits of left- and right-hand 
sides of the equality (14.611 as x approaches oo. The interval [l + ^,2 + i] contains 
a neighborhood of 2 + e, and the continuous P(x) is bounded on the compact set 
[l + |,2 + i]. Indeed (13.81) means ||P|l^oo < for some C > 0. Then, the equality 
(|4.6|) implies 

lim P(x) = 0. 

Next, we consider the positivity of (—l)"'P^")(x) for n = 0,1, 2, 3,4. Similarly with 
the precious lemma, we first study the properties of the reflection only with respect 
to a single inclusion. Second, we apply such properties to the case of two neighboring 
inclusions. 

According to plan, we consider the properties (14.811 . (I4.9|l . (I4.1()L (I4.11|l and (14.1211 
of the reflection with respect to a single inclusion Bq that denotes the unit sphere 
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i?i (0,0,0). As defined before, Ro{h) denotes the reflection of a given harmonic function 
h with respect to Bq. Suppose that for n = 0,1, 2 , 3,4, 

(9"9y/i(a:, 0 , 0 ) > 0 on ( 0 , 1 ]. 

We shall show that for n = 0,1, 2, 3,4, 


{-l)"-d^dyRo{h){x,0,0) > 0 on ( 0 , 1 ]. 

This means (14.8L (14.91) . (I4.10L (14.111) and (14.121) . 

To do so, we use the equality ( 021 ) in the previous lemma or (EH). First, the 
positivity of dyRo{h){x,0,0) results from the equality immediately as 

9.yi?o(h)(a:, 0, 0) = , 0, 0 I — — [ sdyh (s,0,0) ds 

x-^ \x J X Jo 

— ~J s , 0 , 0 ^ — i9y/i (s, 0 , 0 )^ ds > 0 (4.8) 


due to the increasing property of dyh {x, 0, 0) > 0. Second, dealing with the decreasing 
property of dyRo{h){x,0,0), we take a derivative of (14.7|) and then, the increasing 
assumption of dyh{x, 0 , 0 ) yields that 


dxdyRo{h){x,Q,Q) 






sdyh (s, 0 , 0 ) ds — 


0 





dyh (s, 0 , 0 ) 


ds < 0 . 



(4.9) 


Thus, dyRo{h){x,Q,Q) is decreasing. Third, the concavity result can be also obtained 
in the same way. Thus, 


d1dyRo{h){x, 0,0) 

>^Jjs (^Udyh Q,0,0^ - 2dyh{s,0,0)^ ds > 0. (4.10) 

Fourth, we have similarly 

(2^ 

-,0,o]+-^ f sdyh {s, 0,0) ds < 0. (4.11) 

X j Jo 

and 

, 192 /I \ 24 fi 

d^dyRo{h){x,0,0) >—I —, 0,0 j-r / sdyh{s,0,0) ds > 0. (4.12) 

X \x J x° Jo 

At last, we are ready to prove this lemma. By Lemma l3.11 


u(x) = iJ(x) + Ri{H){k) + i? 2 (id)(x) + i? 2 (i?l(id))(x) + i?i(i? 2 (id))(x) + • • • , 


and 


dyuix, 0,0) =dyH{x, 0,0) + dyRi{H){x, 0,0) + dyR 2 {H){x, 0,0) 

+ dyR2iRiiH)){x, 0 , 0 ) + dyRi{R2{H)){x, 0 , 0 ) + • • • , 
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where i?i and i ?2 are the reflections with respect to the insulated inclusions Bi and 
i ?2 as defined in (13.51) . We apply (14.8L (14.911 . (I4.10|) . (14.1111 . (14.1211 . Let n = 0,1, 2, 3,4. 
Since dyH{x,0,0) = 1, 


{-l)^d2dyRi{H){x,0,0)>0 fora;>-| 

and 

d'!^dyR 2 {H){x, 0 , 0 ) >0 for a: < ^. 

In the same way, one can show by the mathematical induction that 


{-ird:idyiiRiR 2 r{H))ix,o,o) > o, 
{-ird:dy{Ri{R 2 Rir{H)){x,o,o) > o 

a:a,((i?2i?ir(i/)))(x,o,o) >0, 

d:idyiR 2 iRiR 2 r{Hmx,o,o) > o 

for X < ^, when m G N. Two representations (EH) and (EH of the solution u yield 
0) =dyRi{H){x,0,0) 

+ dyRi{R2{H)){x, 0 , 0 ) + dyRi[R2{Ri{H))){x, 0 , 0 ) + • • • , 

and by (|3.10ll . 

(-l)”5”5yi?i(iL) (a;,0,0) > 0. 

We thus have 

(-l)”5:aj,i?B,{a:,0,0)>0 

for X > —I that implies this lemma, since P{x + 1 + |) = dyRsi {x, 0, 0). □ 

Another property of P is provided by the following lemma based on the previous 
lemma. This property is used to reduce eh into an ordinary differential equation in 
Lemma 14.91 

Lemma 4.4 For n = 1, 2, 3,4, 


(x-l)”p(”-i)(x) <(x-l)P 


(x—\ 


\ 2 


+ 1 


/or any x > 1 . 

Proof. For n = 1,2, 3,4, the decreasing property of in (l,oo) is provided in 

the previous lemma. In the case of n = 1 , it yields 

|(x - l)P(x)| < (x - 1)P 


for any x > 1 . 

Let n be one of 2,3,4. By the mean value theorem, for any x > 1, there exists 
xo S (x — , x) such that 



on—1 

P<'^-^^(x 1=4) p(”- 2 )(x) 

< 

(x f=4) 


X — 1 

X — 1 
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since the value of (—!)*■" is always positive. It follows from the decreasing 

property of that 




< 


P("-2) {x - 


X — 1 


When n = 3 or 4, we continue this process so that 


p(”-i)(x) 


< 


p{n-2) 

X — I 


pin-3) _ pin-3) 

X — I 

pin-3) _ 1 ^) 

(X- 1)2 

^ P(x-^)-P(x-^) 

^ (x 

(x-l)"-i 


< 


k-il 


Thus, we have this lemma. 


□ 


The fundamental equation (14.61) can be rewritten as 


1 1 
2 ^ 


sP(2 + e — s)ds = xP(x)- -P ( 2 + e-) . 


The left-hand side is positive by the following lemma. The value of the left-hand side 
is very important, since the blow-up rate of P is proportional to 


sP(2 -I- e — s)dsj at X = 1 -I- 2i/e. 

Refer to Lemma [4.15! for the details. 

Lemma 4.5 

sP(2-I-e — s)ds <-— (4.13) 

2 x 2 

for any x € [1, 2 -|- e]. 

Proof. First, we show that the inequality (14.131) is valid for x on the restricted interval 
[ 1 ,P 2 ]- Here, p 2 is the fixed point given in (14.51) . Second, the inequality on [1, 2 -|- e] is 
proved by contradiction. 

According to plan, we prove that 

sP(2 -I- e — s)ds — ^4- < 0 for any x G [ 1 ,P 2 ]- (4-14) 

2 x 2 

It is easy to show that x < 2 -|- e — 4 and x > 1 for any x G [ 1 ,P 2 )- The decreasing 
property of P in Lemma 031 yields 0>P(2-|-e — 4)— P{x) and P (2 -|- e — 4) > 0 
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for any x G [ 1 ,^ 2 )- By Lemma 0^ 


0 >x —P 2 + e- - -P{x) 


1 1 


= / sP(2 + e- s)ds- --y 


for any x £ [ 1 ,^ 2 ], since x > 1 and 2 + e —— = P 2 > 1- Thus, we got the result (14.141) 
restricted on [ 1 ,^ 2 ]- 
Suppose that 

/ sP(2 + e — s)ds — - —2 > 0 
Jo 2 Xq 

for some Xq G (1,2 + e]. By the mean value theorem, there exists a point sq G (0, —) 

1 

such that P(2 + e — So) > 1, since sds = The decreasing property P yields 

P(2 + e - s) > 1 

for any s > sq. For any x G [1, xq], y > ^ > sq so that 


11 

sP{2 + e-s)ds-- — 


sP(2 + e — s)ds — 


1/1 1 


0 V t-2 rf.2 

Z \ X Xg / JQ 


sP(2 + e — s)ds — g —2 
Z Xn 


> [ sds —\\\ - ^1 = 0 . 


O \ rr2 rr.2 

^ \ •^ “^0 


This leads to a contradiction for the first result (14.141) in this proof. Thus, we have 


sP(2 + e — s)ds — 7 :^ < 0 
2 x^ 


for any x G [1,2 + e]. 


□ 


Remark 4.6 It follows from Lemma that 


0 < - — / sP(2 + e — s)ds 

2 ./n 


for any x G [1, 2 + e]. T/iits, we have 


P {1 + s)ds < 1, 


and the decreasing property of P yields 

nl 


J P (1 + s) ds < 3 


due to the smallness of e. 


16 






The function P is a solution to a second-order ordinary differential equation as 
will be seen in Lemma Td.91 The solution P is decomposed into a particular solution 
and a linear combination of two homogenous solutions. The following lemma is used 
to estimate the coefficients in the linear combination that are essential to find the 
blow-up rate of P. Refer to Proposition 14. 131 for the details. 

Lemma 4.7 

371/eP (1 -I- J^/e - 7^e) -I- (^(7^ - l)e - P'(l -I- 7v^) > 0 
*”2<7< To^- 

Proof. Applying Lemma [4.51 to dSl), we have 

0<P(x)-±p(2 + .-i) 

We estimate the ingredients in (14.151) . Let x = H-y-^/e, while 2 < 7 < Then, 

1 - ^ < 3-fy/e, 

since 7 < Note that a; > 2 -|- e — i, since x > p 2 due to the conditions 7 > 2 

and (1431) . The mean value theorem provides the existence of xq G (2 -|- e — y, x) such 
that 


P(x) — P(2-|-e-) = (x-l -2 — e) P'(xo) 


< (x-l -2 — e) P'(x) < 0, 


since P'(xo) < P'(x) < 0 by the monotonic property of P' in Lemma 021 Since 
2 < 7 < Yg^, we also have 

- (^x+^-2-e^ <- (( 7 ^ - l)e - 7 ^e^) < 0 

and 2 -l-e—i>l-|- ^^/e — 7 ^e that implies 

0 < P (^2 -H e - < P (1 -H 7 \/i - 7 ^e) 

due to the decreasing property of P. Applying these bounds above to (j4.15|l . we have 
this lemma. □ 


We consider the property of 2 -|- e — ^ in the equation (14.61) , since the equation is 
the key ingredient in the proof of the first main result. 

Lemma 4.8 Suppose that the sequence {x„} satisfies 

Xi = 1 
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and 


Xn+i = 2 + e- for n G N. 

Xn 


Then, for any n < 
and 

Proof. One can show that 

where 

as given in (Ii3|) , 


= 1 + (« - l)e + o„ 

|o„| < lOn^e-v/e. 


Xn=P2 + 


2 + e — ‘ 2 p 2 
cod"-i + 1 ’ 


Co = 


P 2 = 1 + y e + 

1 - P 2 + e 


©■ 


+ 


and d = 


2 

P2 


1 -P 2 

We estimate p 2 , cq and d in (14.161) . Thus 

1 1 


2 +e-p2 


P2 = 


e+-e+-eV^+0{e^+i), 


1 


Co — 1 — \/e + 2^ > 1 ~ a/c + 2^ 


and 


(4.16) 


d — 1 + 2\/e + 2e + (^^(ev^) > 1 T 2\fk + 2e 
for small e > 0. It has been proved in that 

1 + (n — l)a; < (1 + x)^~^ < 1 + (n — l)a; + (n — l)^a;^, 

supposed that x G (0, 2) and (1 + < 2. Since n < cod"”^ can be estimated 


as 


1 + (2n — 3)\/e < cod" ^ < 1 + (2n — 3)\/e + lOn' 


e. 


Applying these bounds above to (I4.16L we have 


1 


-e—•v/e(l+ - — n \/e 


where 


= 1 + (n — l)e + o„ 

|o„| < lOn^e-v/e. 


□ 
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4.2 Proof of Propostion 14.11 

The decreasing property of dyRsi (x, 0, 0) was presented in Lemma 1431 We first show 
the existence of a constant rg regardless of e such that 

P(l + ro^)<^^. 

e 2 


Second, a relation between the values P(l) and P{1 + rgi/e) is established in (14.4011 . 
Then, we can prove that 


m< 


1 


2-V2 

e 2 


This implies the first bound (14.11) in Proposition 14.11 The second bound (14.21) is 
also presented by virtue of the positivity and decreasing property of dyRsi (x, 0, 0) in 
Lemma 1131 


4.2.1 Estimate for P(1 + rQy/e) for a large tq > 0 

For the sake of convenience, we begin by defining t and / as 

J t = X — 1, 

\ f{t) = Pil + t) = Pix) 


for t>0. 

Speaking of the scheme, the function /, defined in [0, 00 ), is a solution to the 
ordinary differential equation in Lemma 14.91 The function / can be decomposed into 
three functions in (I4.21|) as follows: 


/ — /p + Cafa + Cpfp, 


where fp is a particular solution, and fa and fis are two homogeneous solutions sat¬ 
isfying 

for t > lOy/e. The boundedness of fp is provided in Lemma [4. 101 and the boundedness 
of Ca and the smallness of Cp are also derived by Lemma r4.11l Hence, we can estimate 

P(1 -b roVe) = f{roy/l) < fairoV^) - 

e 2 

in Proposition 14. 131 and Remark 14.141 This is the scheme to estimate for P(1 +rQy/e) 
for a large rg > 0. 

In the following lemma, we establish the ordinary differential equation which / 
satisfies. 

Lemma 4.9 

“ e)/"W + 5</'(t) -b “2fit) = + ffW 


and 


for any t > lOy^. 


I5WI 


< 


tf[\ 
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Proof. By Lemma [4.21 


11 1 / 1\ /■- 

=-yP I 2 + e-) + / sP{2 + e - s)ds + xP{x) 

2x^ x^ \ xj Jq 

for any a; > 1. Taking derivative, we have 

- ^P' f2 + e - -) + xP'{x) + P{x) + f2 + e - 
X'^ \ X J X'^ \ X J 

= - ? - i) - + (-?+") ^' <"> 

1 + P{x) + -4f^’f2 + e-i)- P(a;)) . 


(4.18) 


Since 2 + e—i=a;+^2 + e— 
mean value theorem yields 

(t2-e)/"(t)+5t/'(t) + 2/(t) 


= X + — + 0{t^) in 0 < t < 1, the 


1 


(l + t)3 


+ (3t + + (—1 + 3< + 0{t^)){e — + 0(t^))f'{ti) 

h) + 4t(e + 0{t^))f"it) + O ((e + ff) /'"(h) 

+ o{t)f{t) + o{f)f{h) + o{f)f{t 2 ) + o{t^)f"it) + o{t^)rit3) 


(4.19) 


for t > %/k, where the points ti, ^ 2 , ^3 are located between — and t, and they 
are depending on t. Lemma 14.41 means that 


for t > 0. Thus, 


for any t G (lO-v/e, i). 


t"/(”-i)(t) 

mi < 




tfil 


□ 


Now, we consider the solution to 


{t^-e)r{t) + 5tf{t) + 2 f{t) = -- 


1 


(l + t)3 

We shall find three proper functions fp, fa and fp that satisfy 

1 


+ g{t) for t > IQy/l. 


and 


- ^)fp (t) + ^tfp{t) + ^^3 + 9{t) 


{t^-e)f''{t)+5tm+mt) = o 


(4.20) 


for i = /3. The general solution is decomposed into the three functions as follows: 

/ = /p + Cafa + Cpfi^, (4-21) 
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where fa and fp is homogeneous solutions defined as 


oo n 


n=l 


n—1 k—1 


(2A: - y2)(2fc + 1 - v/2) 

2k(2k - 2 V 2 ) 

" A (2fc + V2)(2fc + 1 + ^2) 


2 fc( 2 fc + 2 y 2 ) 


(4.22) 

(4.23) 


for t > lO-v/e. 

The functions fa and fp can be established by induction. To do so, we regard 
fa and fp as the sums fl^n , where fao = and fpo = 

^- 2 -V 2 solutions to t'^fioit) + + 2 /io (0 = and /„ is the solution to 

+ 5i//„(i) + 2/jn(t) = for i = a, and n = 1, 2, 3, • • •. The functions 

fa and fp are dehned well on [lO-y/e, 00 ), because 


fanjt) 
fa{n—l) (f) 

fpnif) 

fP{n—1) {t) 



(2n - V2)(2n + 1 - \/2) 


2n{2n — 2\/2) 


{2n + V2){2n + l + V2) 

VI 

U 2 J 

2n{2n + 2y/2) 


and the variable t > 10^/e. Moreover, we have 


fa ~ t and fp^t 


- 2 -V 2 


Dealing with (|4.21|) . we consider the contribution of fp to /. The boundedness of 
fp is derived in the following lemma. 

Lemma 4.10 There are a particular solution fp to (I4.20p and a constant Cq such 
that 

I/pWI < 1 

and 

\fp{t)\ < \ 

for any t > Coy/e. 

Proof. We shall find the sequence of functions {fpn} satisfying 

ff^oit) + 5t/;o(i) + 2fpoit) = G{t), (4.24) 

+5t/;„(t)+2/p„(t) = e/;(„_i)(t) (4.25) 

for n = 1, 2, 3, • • •, where 

The sum /p” desirable function fp. Without any loss of generality, we 

assume in this proof that 

\9{t)\ < tf (4.26) 

for any t > 10\/e. This is a simplication of the inequality | 5 (t)| ^ tf (^ in Lemma 
HU for convenience. 
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Define fpo{t) as 


fpo{t) = 


1 


f-2+V2 


lOy/e W 


1 -2^2 


'lOV^ 


G{s)s^-'^dsdw. 


Then, fpo(t) is a solution to (14.241) . Moreover, we estimate f^{t) 
Remark 1121 / 2 e/ (f) ^ 2 . Lemma and (14.261) yield 


(4.27) 
for n = 0,1, 2. By 


\fpo{t)\ < 


f2+V2 


D-V 2 . 


'lOV^ W 


l-2\/2 


10^ 


/(i) 


, 2 -A 


dsdw 


1 


< ”1 -F / ^ 

- 2 t2+V2 


l + \/2 


'lov^ 


/(I) 


dsdw < 2 


in lO-y/e < t < 1. Taking the derivative of (14.271) . we can get similarly 
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and by ()4.24l) 


in 10^/e < t < 1. 

We also define fpn{t) as 

fpn(t) = 


\fpo{t)\ < - 


|/;'oWI<^(60+|5(i)l) 


1 


^ 2 +V 2 


/l^ W 


I- 2 V 2 


/ 120 e 


,^fp{n-l)^ 


I-V 2 


dsdw. 


(4.28) 


Then, /p„ (t) is the solution to (14.251) . In the same way, we can prove by mathematical 
induction and (|4.26p that for any n = 1, 2, • • •, 


\fpn{t)\ < |/;„(t)| < ~ (4.29) 

and 

\^fpn{t)\ < ^^(60 + mi) < ^(1 + l5(t)l) (4.30) 

that is the right-hand side of (j4.25l) . while \/120e < t < ^. Hence, the sum /po 

is well defined and is the desirable function fp. □ 


We shall consider the contributions of fa and fp to / in Proposition 14.131 since 
the boundedness of fp was derived in the previous lemma. To do so, we need Lemmas 
I4.11l and and 14.121 

Lemma 4.11 

f (7\/e - 27 ^e) + (( 7 ^ - l)e - 7^e^) > 0 

*« 2 < 7 < 15 ^. 

The lemma above is a rewritten version of Lemma 14.71 since f(t) = P{1 + t). 

We also need the following lemma to prove Proposition 14.181 
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Lemma 4.12 Suppose that the constants M > 0, Ca, C /3 > 0 and Cq > 0 satisfy 



M + Ca ^ ^ : 

t2-V2 - 

Then, 


Proof 

For any t > Coy/^, 


Mt^- 


- f2+V2 ^ - ^^0^- 


Let t = s{Co^/e). For any s > 22 v^, Then, 


Mt^-'^ + Co, > M{ 


\ 2 -V 2 ~ ~ ~ 

(coVi) + a > Cp— —^ > 2Cp 


(c'ov^) 


2\/2 “ ^ ^ 2\/2 ■ 


□ 

Using Lemmas I4.11l and l4.121 we estimate fit) in [ 7 o\/e, in the following propo¬ 
sition, supposed that 70 is sufficiently large regardless of e. We shall prove in Subsec¬ 
tion HTHIll that P (1 -I- 7 o\/e) = / (7o-\/e) has the same blow-up rate as P(l). In this 
respect, the estimate for / ( 7 o\/e) is meaningful. 

Proposition 4.13 There is a constant 70 independent of e such that 

fit) < Ut) ^ 


for t e [ 7 ov/e, ^). 

Proof. We consider the constants Ca and Cp in the decomposition of / as 


/ — /p + Cafa + C/sf/S 


(4.31) 


as given in (14.211) . We recall the behavior of fp, fa and f/ 3 . The boundedness of 
I1 /p(OIIloo([CoU£ i]) ll^•^p(^)|lL°o([c Ve -]) presented in LemmaSTTUl From the 

definitions (14.221) and (14.231) . 


/-(*) = 7 ^ (1 + o (^)) “^ ^ (1 + O (i)) 


(4.32) 


for t > Coy/e, supposed that Cq is sufficiently large regardless of e. In this proof, we 
consider two cases when < 0 and C /3 >0, respectively. 

In the first case when C /3 < 0, the decomposition (14.3111 and the positivity of / 
in Lemma 1221 yield that fp + Cafa > —C/ 3 f /3 > 0. By the boundedness of fp, there 
exists Ml > 0 such that 


2-^1 > ll/p(i)llL=o('[c.o^,i]) ■ 


(4.33) 


Then, 


Ml + Cafait) > -Cf3f/3{t) 
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for t > Co-v/e. Note that fa — ^ 2 ~V 2 fp — ^ 2+^2 ■ Lemma 14.121 there is ro > 0 
such that 

i (Ml + Cafait)) > -Cfifp{t) (4.34) 

for t > ro-v/e > Coy/e. Thus, 

1 ^ f + Midt 

J To y/l 

> J -i^Cafa + ^ 2 “^^ (^2 

n 1 

> C,. / -fadt > Ca, 

Jro^/f. ^ 

1 

since 1 > Jf^^fdt by Remark |T6l Hence, we use (14.341) to get 

fit) <1 +Cafait) < fait) 


for t > roy/e. 

In the second case when C /3 > 0, we have f > fp + Cafa- Then, it follows from 
Remark 14.61 and Lemma [4.101 that 

Ca [ fadt < f f - fpdt < 1 . 

J Cq^/e J Cq^/e 

Thus, since Jc^^fa^t — 1, the constant 

Ca < 1 . 

On the other hand, we estimate Cp. Note that 3 — (2 — y/2) > 0 and 3 — (2 + ^/2) < 
0. Applying (14.311) and (14.321) to Lemma 14.111 we can find a positive constant M 2 
regardless of e so that 

+ (4.35) 

and 

M 2 + Cafait) >Cpfpit) (4.36) 

in Co\/e < t < jh, where the boundedness of fp is also used to get M 2 . Thus, we have 

0</(t)<M2 + C„/„(t)</„(t) 

for t G [Coy/e,^]. 

Hence, regardless of whether Cp > 0 01 Cp < 0, 


in roy/e <t<^. 


fit) < fait) 


□ 


Remark 4.14 It follows from Proposition\4-13\ that 


Pi^ + t) < 


1 


in rgy/e < t < ^. Hence, we have 

P(l + rov^) < 

e 2 
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4.2.2 Estimate for P{1) 

The estimate for P(l) results from P(1 + which was presented in 

Proposition 14.131 and Remark 14.141 We first prove Lemmas 14.151 and 14.161 to show 
P P TMV^) ^ 2 ^V 2 second establish the chain (14.401) of correlations between 

e 2 

P(l) and P (l + Y^-\/e)- Hence, we can obtain the estimate for P(l) as follows: 

^( 1 ) < 


2-V2 ' 

e 2 


Lemma 4.15 


0 < - - 


1 1 

2 - / {2 + e-t)P{t)dt<^. 


(4.37) 


Proof. Let x, = 1 + 7 o-v/e where 70 was given in Proposition 14.131 and rp > 3. 
Since cc* > 2 + e — the decreasing property of P in Lemma [4.31 implies P(xP} < 


P (2 + e — By Lemmas 14.51 and 


0 < 


1 1 
2 ^ 


sP(2 + e — s)ds = X, ( P(x*)-^P ( 2 + e- 


By Proposition 14. 131 and Remark 14.141 

0 <t Pi^xP) — —tP I 2 + e — — 


-P(x,)-p(2 + e-^) + (l-^)p(^2 + e-^ 


J \ X 

By Remark 14.61 and the positivity of P in Lemma na 

f2+e 

'1+2^ 
f 2+e 

J 2 +e-^ 

since 1 + 2^/e < 2 + e — -P due to xq > 3. 

Lemma 4.16 


2-\/2 1-^2 ’ 
6 2 6 2 


1 f 

0<-— / (2 + e — t)P{t)dt 

2 Jl+2^ 

^ 2 +£ 2 
/ {2Pe-t)P{t)dt + rayfe<^^, 

d 2 +e- 4 - f 2 


1 1 

- 2 ^ 


2-^2 


6 2 


Proof Applying Xa = 1 + T^\/e to Lemma IL^ we have 


□ 


(4.38) 


i — xP [ sP{2 + e — s)ds = xpP[x) — P [2 + e - 


and note that 


Xa 7 2 + e — — < 1 + 2\fk. 
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It follows from (14.3711 that 


^ ^ / “ sP(2 + e - s)ds 

do 

1 


(4.39) 


= Xa^P{xa) - P 1^2 + e- 

> (Xa^ - l)P(Xa). 


Then, 




2-V2 

e 2 


□ 


To estimate P(l), we establish a correlation between P(l) and P (l + Y^-\/e)- 
The correlation is a long chain based on 


P ( 2 + e-) = x^P{x) +x^ sP{2 + e — s)ds — - 


in Lemma 1321 In Lemma IT21 the sequence {xn} was defined as 


Xi = 1 and Xn+i = 2 + e-for n G N. 

Xr), 


Then, 


P(x„+i) = xlP{x„) +; 


sP(2 + e — s)ds — 


for n G N. Let uq = 


20^/7 


. Then, 


no —1 


P {Xno + l) + ^ ^ ^n+1 ■ ■ ■ ■ (Xn+l) 


n—1 


~ ^no'I^(^no) + • 


sP(2 + e — s)ds — 


no-l / 

+ ^n+l xlP{Xn) + : 

n=l \ 


sP(2 + e — s)ds — 


By cancellation, 


2^1 • •••• a^no 1^(1) 


= P(x„„ + i)+ I 2 “^no 


sP(2 + e — s)ds 


no —1 


^n+l ' 


— x^ sP{2 + e — s)ds . 


(4.40) 


It follows from Lemma [4.81 that I + < Xn^+i- Thus, by Lemma [4.161 and the 

decreasing property of P in Lemma 14.31 


p (a^no+i) < P + ^ 


2-n/2 

e 2 
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We use Lemma 14.81 again so that 

and ccra > 1 for n € N. Note that 1 + 2i/e > p 2 > for n = 1, • • • , ug, where p 2 is 
the fixed point in (14.51) . By Lemma 14.371 

1 1 o 

> 7 ;-Xn / sP{2 + e - s)ds 

for n = I,-- - ,no. Therefore, the bound (14.401) can be reduced into the desirable 
bound (14.11) as follows: 

dyRB, (-|.0>0)=P(1)<^. 

6 2 

As mentioned earlier in Subsection HU the second bound (03 is also obtained 
immediately by the positivity and decreasing property of dyRsi (x, 0,0) in Lemma 

[121 □ 

5 Proof of Theorem 12.2 

We establish the lower bound in Theorem o under the assumption that 

H{x,y,z) = y in 

The upper bound of Vu and the properties presented in Section m are used to derive 
the lower bound. The directional derivative dyU is decomposed as 

dyu{x, 0,0) = 1 + dyRsi {x, 0,0) + dyRB 2 {x, 0,0) 

for I a; I < |. As defined in Subsection 14.21 

f^t)=P{l+t)=dyRB, 

for any t > 0. 

In this proof, we shall prove the existence of tg > 0 with 

fito) > (5.1) 

e 2 

By Lemma 1121 (x,0, 0) = dyRBi{—x,0,0) > 0 for x < | and (x,0, 0) is 

decreasing for x > — |. Then, dUD implies 

dyu{xo,0,0) > 

e 2 

for some xg € (—f, §)• This is the desirable lower bound. 

To show the existence of to, we need a negative particular solution fp to the 
equation (14.201) . 

Lemma 5.1 Let g(t) and G(t) be as given in (14.201) and (14.241) . There exist a negative 
function fp, three positive constants M, sg > 100 and Sq < ^ such that 

+ ‘2'fp{t) = - ^^3 + git) = G{t) (5.2) 
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and 


-M < fp{t) < 


1 


for anyt G {soy/e,So). Here, M, sq and So are the constants regardless of small e > 0. 


It is worth mentioning that the function fp in this lemma is slightly different from the 
one in Lemma 14.101 


Proof. By Lemma 03] and Proposition 14.131 there are constants si and such that 

for any t G {siy/e,Si) and Si < 

As defined in Lemma [4.101 we also consider 


fpit) — fpn{t). 


n=0 


In the similar way to (14.2711 and (14.2811 . we define fpn{t) as 


and 


/po(i) = 


fpn{t) — 


f2+V2 ^yl-2\/2 


G(s)s 


I-V 2 


dsdw 


1 


1 


^2+V2 


siy'e W 


I- 2 V 2 


-fPin-1) '^dsdw 


si-v/e 


for any t G (sii/e, ^i) and n = 1, 2, 3, • • •. For any t G (2sii/€, Si), 

1 


—Ml < fpo{t) < — 
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since G{t) < —j- Here, Mi is the constant regardless of small e. 

We can use the mathematical induction to get the analogues of (14.2911 and (|4.3011 . 
By virtue of in (14.3011 . there are constants sq > 100 and S'o < such that 


\fpn{t)\+t\f'{t)\+t^\f;jt)\< 


1 1 


for any t G (so\/e, Sq) C (si\/e, aS”!) and n = 1, 2,3, • • •. Hence, fp is the solution to 
(EH) and satisfies 

-M < fpit) < - — 

on the open interval (so\/e. So). Here, M is the constant regardless of small e. 

□ 


Lemma 5.2 Let sq o.'^d So be as given in the previous lemma. Then, 


pSo 

/ f{t)dt > 1 . 

Sq \/e 


Proof. Lemma 14.151 yields 


1 

'7 Si (2 + € — t)Pit)dt 

4 Jl+2^ 
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for small e > 0, since P{x) < Thus, 

£ 2 



for small e > 0. The decreasing property of / in Lemma 14.31 yields 


10 y 1 + e — SQy/e j 


< 





Thus, we have this lemma. 

Now, we recall the decomposition (14.211) as 

f = fp + Cafa + C'^//3 


□ 


(5.3) 


where fa and fp are defined as (I4.22|) and (14.231) . We consider two cases when Cp <0 
and when Cp > 0, separately. 

In the first case when Cp < 0, the negativity of fp in Lemma |5.II and Lemma |5.21 
yield 

Sq Sq Sq 

Ca j fadt > f Cafa + (fp + Cpfp) dt = f fdt>l 

J SQy/e J SQ\/e J SQy/e 

due to the positivity of fp in (14.231) and (I4.32p . Here, sq and are the constants in 
Lemma 15.11 From the definition (14.221) . J^°^fadt < 1 so that 

Ca > 1. (5.4) 

On the other hand, the positivity of / in Lemma l4.3l and the negativity of fp in Lemma 
EUlyield Cafa{t) > -Cpfp{t)-fp{t) > -Cpfp{t) > 0 on the interval (soy/e, So). Then, 


> 1 , 


Cafa{2so^/e) 

Cpfp{2so\^) 

and considering the defintions (14.221) and (I4.23L 

since the constant sq > 100. Thus, we use the definitions of fa and fp again to get 
Cgfait] ^ i '^ > 2 (5 5) 

for t > Ssoi/e. Applying (15.51) to (15.31) . the inequality (15. 4|) implies 

f(8soVe) > ^Cafa(8soVe) -M> 

when Cp < 0. Here, the constant M is given in Lemma ISTTl 
In the second case when Cp >Q, we shall prove that 


Ca > 1. 


(5.6) 
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By (14.351) and (14.3611 in the proof of Proposition 07131 there are positive constants Ca, 
S 2 and S 2 such that 


Ca M +Co 


1 


t 2-\/2 


>Cp 


1 


t2 + y/2 


(5.7) 


and 


CA{M + CMt))>C^flsit) 


for any t G {s 2 \/e, S 2 ) C (so\/e, 5'o), while the constant M was given in Lemma [5.11 
To prove (15.611 . we subdivide this second case into two subcases, due to (15.711 . First, if 


CaM > -Cp- 


1 


2 ^ {2s2yfeY+'^' 

then, the definition of fp in (14.231) yields 

f‘5'2 

^ 2 s 2 \/c 

In the same way as Lemma 15.21 the negativity of fp in Lemma 15.11 and (15.31) yield 

rS2 rS2 

' 2 s 2 \/e 
r5'2 


Cp fpdt < 

J2so\/e 


1 < / fdt<Co [ ^ Udt + 

J 2S2.yJ~^ 2S2\/~^ 


Thus, since / 2 s 2 ^/a'^t ^ 1) we have 


> 1 . 


Second, otherwise, if 


CaM < -Cp- 


2 ''( 2 s 2 v^) 2 +^’ 

then ([53 implies > Cp For any t > 2s2\/e 


Co-^>Cp- ^ 


' I 2 -V 2 ^ f2 + y/2 ' 


Thus, 


S 2 S 2 S 2 

Ca f fadt > f Cafa + Cpfpdt > f fdt>l 

J 2s2\/^ «/2s2\/c «/2s9\/c 


^82 
/ 2 s 2 \/e 
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2 S 2 


since fp < — 3 ^- This means 
Hence, it follows from (15.311 that 


Ca > 1 . 


/(2s2V^) > Cafa{2s2V^) -M> 

when Cp > 0. 

Therefore, as a result, there is a point to with 


2 -V 2 

e 2 


/(to) 


> 


1 


2-^/2 


e 2 


regardless of whether Cp > 0 or Cp < 0. This is the desirable (15.111 . 


□ 
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